ROME1-1222/9J 



Heavy-heavy form factors and generalized factorization 



M. Ciuchini 1 , R. Contino 2 , E. Franco 2 , G. Martinelli 2 
1 Dipartimento di Fisica, Universita di Roma Tre and INFN 
Sezione di Roma III, Via della Vasca Navale 84, 1-00146 Rome, Italy. 
2 Dipartimento di Fisica, Universita "La Sapienza" and INFN, 
Sezione di Roma, P.le A. Moro 2, 1-00185 Rome, Italy. 



Abstract 

We reanalyze B — > Dtt and B — > KJ/ijj data to extract a set of parameters which 
give the relevant hadronic matrix elements in terms of factorized amplitudes. Various 
sources of theoretical uncertainties are studied, in particular those depending on the 
model adopted for the form factors. We find that the fit to the B — > Dtt branching 
ratios substantially depends on the model describing the Isgur-Wise function and on the 
value of its slope. This dependence can be reduced by substituting the BR{B — > Dtt) 
with suitable ratios of non-leptonic to differential semileptonic BRs. In this way, we 
obtain a model-independent determination of these parameters. Using these results, 
the B — > D form factors at q 2 = M 2 can be extracted from a fit of the BR(B — > Dtt). 
The comparison between the form factors obtained in this way and the corresponding 
measurements in semileptonic decays can be used as a test of (generalized) factorization 
free from the uncertainties due to heavy-heavy form factor modeling. Finally, we 
present predictions for yet-unmeasured Dtt and DK branching ratios and extract fo s 
and f D * from B DD S decays. We find f Da = 270 ± 45 MeV and f D , = 260 ± 40 
MeV, in good agreement with recent measurements and lattice calculations. 
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1 Introduction 



A problem of utmost importance in B phenomenology is the computation of the 
hadronic amplitudes: in recent years it has been realized that the full determination 
of the unitarity triangle from B decays can hardly be carried out without an accurate 
knowledge of these quantities ^] . Unfortunately, the computation of hadronic am- 
plitudes requires an understanding of low-energy strong interactions which is missing 
at present. Even a non-perturbative approach based on first principles, such as lattice 
QCD, fails in computing decay amplitudes involving two or more hadrons in the final 
state [|. 

In the absence of rigorous methods, some simplifying approaches have been de- 
veloped. The most popular one consists in the factorization of matrix elements of 
four-fermion operators in terms of local products of two currents. In this approach, 
the original matrix element is computed as product of the matrix elements of the two 
currents. Attempts to give theoretical soundness to this procedure in the framework 
of the 1/N expansion and of the Large Energy Effective Theory (LEET) can be found 
in refs. [Q, ||. Unfortunately, there are many problems in both these approaches and 
their applicability to exclusive decays is questionable. Independently of any theoretical 
prejudice, there is a priori no reason for this approximation to be accurate in the case of 
B decays. Indeed, none of the expansions developed so far was able to compute correc- 
tions to the lowest order results and estimate the size of the errors. On the other hand, 
the importance of controlling the theoretical uncertainties calls for some phenomeno- 
logical approach to test predictions obtained using factorized amplitudes. To this end, 
a popular method consists in reducing the Wick contractions of matrix elements to few 
topologies using Fierz transformations and color rearrangement. Then, the remaining 
amplitudes are factorized and expressed in terms of the appropriate decay constants 
and/or form factors. In this procedure, some phenomenological parameters are intro- 
duced in order to account for possible deviations from factorization 0,0,0. These 
factorization parameters, denoted as FP in the following, are meant to be fitted to 
experimental data. 

In this paper, we introduce a parameterization of the hadronic matrix elements that 
extends the one of ref. Jl| and allows the computation of the hadronic amplitudes rele- 
vant to Cabibbo-allowed non-leptonic B decays in terms of factorized matrix elements 
and of three real FP. We find that, in the fit of the BR(B -► Dtt) and BR(B -► KJ/V), 
there is a strong interplay between the values of the FP and the model used for the 
heavy-heavy form factors, more specifically on the Isgur-Wise (IW) function and its 
slope p 2 a. This implies that factorization tests are obscured by our ignorance on the 
values of the form factors in the kinematical region relevant in non-leptonic decays 
(q 2 <C qmax)- The model dependence is drastically reduced by using, in the fit, suitable 
ratios of semileptonic and non-leptonic BRs (to be introduced below) instead of the 
BR(B — > Dtt) alone. In this way, we are able to extract (almost) model independent 
FP. With these FP at hand, we use then the BR{B —* Dtt) to determine, within a 
given model, the value of p 2 which may be compared, as a test of factorization, to the 
one measured in semileptonic decays. The value of p 2 extracted from the fit depends, 

1 Here and in the following, unless stated explicitly otherwise , B — > Dtt denotes generically a full set of 
decays of B mesons into a charmed and a light meson, i.e. Bd — > tt + D~ , Bd — ► p + D~ , B + — > ir + D*°, etc. 



1 



Form Factor 


LINSR 


NRSX 


MM-) 

V(M 2 ) 
A (M 2 ) 

M{Ml) 


0.56 [0.49-0.63] 
0.68 [0.61-0.73] 
0.59 [0.54-0.64] 
0.59 [0.53-0.64] 
0.59 [0.53-0.64] 


0.57 [0.52-0.63] 
0.75 [0.69-0.81] 
0.58 [0.54-0.63] 
0.56 [0.52-0.61] 
0.54 [0.49-0.58] 



Table 1: Values of B — > D form factors determined by fitting the B — > Dix data with two 
different models, as explained in the text below. The ranges in square brackets correspond to 
variations of x 2 /dof up three times larger than its minimum. We find that form factors at 
q 2 = M 2 already verify the kinematical relation /o = /+ and A% = Aq valid at q 2 = 0. 



however, on the model used for the form factors. Different values of p 2 compensate, 
indeed, for the different dependence of the theoretical form factors on the momentum 
transfer, thus giving the same values for the matrix elements of the weak currents at 
low q 2 . We conclude that the quantities to be compared with the corresponding ones 
in semileptonic decays are the form factors themselves in the region of q 2 relevant in 
non-leptonic decays (e.g. q 2 = M 2 ~ for B — > Dir decays). This is a real test of fac- 
torization, free from model uncertainties. The values of the form factors extracted from 
our analysis are given in table [l]. In principle, one may extract the five form factors of 
tab. [l] independently. However, in our analysis, all the form factors are related to the 
IW function through the heavy quark symmetry. Consequently, the only B — > D form 
factor measured so far at small q 2 , /+(0), already allows a full test of our approach. 
Its experimental value, /+(0) = 0.66 ± 0.06 ± 0.04 B is in good agreement with our 
findings. Measurements of the other form factors entering B — > D* semileptonic decays 
would check the relations enforced by the heavy quark symmetry. 

The determination of the FP also allows us to predict several BRs, including B 
DK, which have not been measured yet. Our predictions are presented in table 

Finally, using ratios of non-leptonic BRs involving D^Ds final states and the FP 
and form factors from the previous fits, we extract the meson decay constants fz> s and 
f D * , obtaining 

f Ds = 270 ± 45 MeV, f D * = 260 ± 40 MeV, (1) 

in good agreement with recent experimental measurements fu g = 250 ± 30 MeV || 
and lattice determinations, f Ds = 218±f2 MeV (quenched), f Ds = 235tf 5 1ll 7 MeV 
(unquenched) Jul and fo* = 240 ±20 MeV (preliminary quenched) pl| . We also study 
the contribution of charming penguins JlJ and discuss their effects on the predictions 
for the decay constants, which we find non-negligible. 

The paper is organized as follows. In sect. ^ we introduce our FP and two different 
models for the form factors to be used in the phenomenological analysis. Section || 
contains the main results of our fits to the B — > Dtt and B — > KJ/^> branching 
ratios, namely the determination of the FP, the analysis of their p 2 dependence and 
the extraction of the B — > D and B —* D* form factors at q 2 = M 2 . The results of 
these fits have been used for the predictions of yet-unmeasured BRs, including many 
B — ► DK modes. Finally, in sect. ||, we analyze the B — > DD S modes and extract fjj a 
and fr>* , giving an estimate of the theoretical error which includes charming-penguin 
effects. 
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Channel 




LINSR 




NRSX 


Experiment 






[BR x 10 5 ] 


[BR x 10 5 ] 


[BR x 10 5 ] 


B d - 


-> 7r u D u 


14 


[1-58] 


10 


[2-27] 


< 12 


B d - 


-> 7T U L>* U 


15 


[1-63] 


13 


[2-35] 


< 44 


B d - 




6 


[1-26] 


7 


[1-19] 


< 39 


B d - 


n j.n 

-> 


17 


[2-71] 


14 


[2-38] 


< 56 


B d - 


-> K + D~ 


22 


[13-38] 


23 


[16-32] 


— 


B d - 


-» K + D*~ 


22 


[14-37] 


22 


[16-30] 


— 


B d - 


-» K*+D- 


53 


[32-90] 


53 


[38-75] 


— 


B d ~ 


-> K* + D*~ 


67 


[43-110] 


64 


[46-88] 


— 


B+ - 


-> K+D° 


35 


[12-54] 


35 


[18-45] 


29 ± 10 


B + - 


K + D*° 


36 


[13-54] 


35 


[17-44] 


— 


B+ - 


-> K*+D° 


67 


[34-100] 


69 


[42-88] 




B+ - 




87 


[46-126] 


83 


[51-104] 




B d ~ 




1.4 


[0.2-5.7] 


1.0 


[0.2-2.7] 




B d ~ 




1.5 


[0.2-6.2] 


1.3 


[0.2-3.4] 




B d ~ 


-» ET*°Z) 


0.6 


[0.1-2.7] 


0.7 


[0.1-1.8] 




B d ~ 




1.7 


[0.2-7.2] 


1.4 


[0.2-3.6] 









Table 2: Predictions of yet-unmeasured branching ratios. 



2 Factorization, FP and form-factor models 

In this section, we present our parameterization of the hadronic amplitudes and discuss 
its relation with other popular choices. We also introduce two different models for the 
form factors used in our phenomenological analysis. 

Consider the matrix element of some composite operator appearing in the AB = 1 
weak Hamiltonian, between the B meson and two final pseudoscalar or vector mesons. 
In general this operator can be written as the product of two currents. If one of the 
currents has the correct quantum numbers to create one of the final state mesons 
from the vacuum, then the matrix element can be factorized. The physical idea is 
the following: the quark pair produced by this current acts as a color dipole, weakly 
interacting with the surrounding color field. If the transferred energy is large, the quark 
pair has no time to interact before hadronizing far from the interaction point |i~2|| . 

As an example, we discuss the factorization of the amplitudes entering the decay 
Bd — > D~ir + . In this case, the two relevant matrix elements (a and (5 are color indices) 

(ZTTT+lQilBd) = {D-ir+\b a ^(l - 75 ) C/3 u p ^(l - j 5 )d a \B d ), 
(D-ir + \Q 2 \B d ) = (D-7r + \b allx (l- l5 )c a u^(l--/ 5 )dp\B d ), (2) 

can be Wick-contracted according to two different topologies, that are usually denoted 
as connected (CE) and disconnected (DE) emissions, respectively. Color indices can 
be rearranged using the algebraic relation, 

botfi&pa = J-;5aa5pi3 + 2^^, (3) 
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where ./V is the number of colors, 5 is the Kronecker symbol and the t a are the SU(N) 
color matrices in the fundamental representation, normalized as tr(t a t b ) = 5 ab /2. Using 
this relation, one obtains 




CE DE octet terms 

In the factorization limit, no gluon exchange occurs between the quark pair of the 
emitted meson and the other quarks, so that the octet terms vanish and the relation 
between DE and CE becomes simply CE = DE/N. Exact factorization is known to 
fail, however, in reproducing D phenomenology |D|. For this reason, it is customary 
to introduce several phenomenological parameters to account for octet terms (and in 
general for the different sources of factorization violation) . These parameters may be 
extracted from the experimental data. An example is provided by the generalized 
factorization of ref. 0. In this case the relevant contractions are rewritten, without 
loss of generality, as 

DE=(l + e 1 )DE fact , CE=(Jj + j^-}(l + ei)DE faet , (4) 

where the two parameters, ei and eg, vanish in the factorization limit. 

In this paper, following ref. |l|], we adopt a different parameterization, given by 

DE = aDE fact , CE = ae i5 ^DE fact , (5) 

where DE and CE are given in terms of three real parameters a, £ and 5^. Note 
that there is no inconsistency between the two parameterizations: in general, there 
are three real parameters, namely two moduli \DE\, \CE\, and one relative phase, 
&rg(CE) — arg(DE). These correspond to our three real parameters a, £ and 5^ or to 
the two complex parameters e\ and es in eq. (Q), one of which can always be chosen 
real. The relation between the two sets of parameters is given by 

«-!+«., e^-^ + T^- <«> 

As recently stressed in ref. [Q, these phenomenological parameters are renormal- 
ization scale and scheme dependent, as much as the original matrix elements, since the 
factorized amplitudes are insensitive to both the scale and the scheme. This depen- 
dence is required to cancel the corresponding dependence in the Wilson coefficients, up 
to the order at which the perturbative calculation is done. Note that, in order to study 
the scale dependence of the parameters, the next-to-leading order (NLO) determina- 
tion of the effective Hamiltonian is required. Being scheme-dependent, any physical 
interpretation of the "factorization scale", namely of the renormalization scale (if it 
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really exists) at which exact factorization holds, is meaningless. Nevertheless, the FP 
can be precisely extracted from data, once the renormalization scale and the scheme 
are fixed. Their values will depend, of course, on these choices. We will use the NDR- 
MS NLO Wilson coefficients computed at // = 5 GeV, as given in ref. ||15|| . In the 
following, it is understood that we determine a, £ and 5^ using this choice of the scale 
and of the renormalization scheme. 

After the introduction of the FP, the only amplitude that remains to be computed, 
namely 

DE fact = (D-n + \hnO- ~ 7s)c vrf(l - -ys)d\B d )\ faet 

= (D-\hn(l ~ j 5 )c\B d )(TT + \u^(l - l5 )d\0) , (7) 

can be easily expressed in terms of the B — » D semileptonic form factors and of the 
decay constant f n . 

In this example, only left-handed currents appear. In general, considering also 
penguin operators, there are diagrams involving different Dirac structures, e.g. 7^(1 — 
75) <8> 7^(1 + 75) and (1 — 75) g) (1 + 75). In the case of interest, the right-handed 
current always appears in the matrix element of the emitted meson, (M\qj^(l+ , f5)q'\0), 
while the current entering the other matrix element is always left-handed. Therefore 
only the vector or the axial current separately contributes, depending on the quantum 
numbers of the emitted meson. Consequently, assuming that both left-left and left- 
right operators can be described with the same set of FP, the relation between the 
corresponding matrix elements becomes trivial. Similarly, matrix elements of operators 
with a (1 — 75) g) (1 + 75) Dirac structure can bfi connected to the current-current ones 
via the vector and axial vector Ward identities a. In summary, using factorization, one 
only needs to compute matrix elements of currents, that can be expressed in terms 
of form factors and/or decay constants. However, these relations among insertions 
of different Dirac structures only hold for factorized amplitudes. By using only one 
set of FP, we implicitly assume that the same relations hold for the original four- 
fermion operator matrix elements. This simplifying assumption allows us to account 
for penguin-operator contributions using factorization, albeit in a model-dependent 
way. 

In our analysis, we use the same FP for i) matrix elements connected by SU (2) fla- 
vor symmetry; it) matrix elements with the same quark content, differing only for the 
angular momentum of the final hadrons. The first assumption has sound phenomeno- 
logical motivations; the second is reasonable since some of the differences among matrix 
elements with pseudoscalar and/or vector meson final states are already accounted for 
by factorized matrix elements. 

For a generic transition B — > P(V) of a B going into a pseudoscalar (vector) meson 
of mass M, momentum p (and polarization e), the form factors in the helicity basis are 
defined as 

(P{P)\%\B{ PB )) = f + (A{p + PB )„ - q , Ml ~ M2 I + /o( g 2 ) g ,- (M| " 



1 q I q 

<y(p,e)\%\B(p B )) = MB 2l +M e^ aP e*»p a p B V(q 2 ), (8) 



2 In this case the amplitudes depend on the quark masses which we take to be defined in the same 
renormalization scheme, and at the same scale /x, as the four-fermion operators. 
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(V{p,e)\A^\B(p B )) = 2MA (q 2 )q^^+A 1 (q 2 )(M B + M)^e;-q^-^ 

; (PB + • 9 - qti^-{M B - M ) 



where and are the vector and axial currents respectively. 

For heavy final mesons, the form factors can be connected to the HQET functions 
see e.g. ref. fig), using the following formulas 



(F(t7,e)|y M |S(« B )) = v/MbM {iM^e^e*^?}, (9) 
{V(v,e)\Ar\B(v B )) = ^Mb~M {uM(y + iK 

-£,A 2 {y){(* ■ v B )v Bfl - U s {y)( e * ■ v b)vh}, 

where v B and v are the 4- velocities of the initial and final meson respectively and 

M % + M 2 -q 2 

V — Vr-v = — . 10 

y 2M B M V ; 

In the heavy quark limit, the functions £i(y) are all related to the IW function 

£+ = Cv = U = U 3 = C(y), e- = U 2 = o, (ii) 

with the normalization of £ fixed by the heavy quark symmetry, £(1) = 1. The 
can be written as 

£<(v) = {«, + + 7<(y)} (12) 

where the functions Pi(y) and 7«(y) take into account the perturbative 0{a s ) correc- 
tions and the 0(l/m) terms respectively B. Following ref. [17|, we used for fh the 
reduced charm-bottom mass i.e. fh = 2.26 GeV. 

We are now ready to introduce the two models that we will use in order to study 
the form-factor dependence of the FP. We denote these models as LINSR and NRSX: 

• LINSR : the first model uses the heavy-heavy form factors defined in eq. (|i~2|), 
taking the (3i(y) from ref. [17] and neglecting the 0{l/m) corrections, i.e. Ji(y) = 
0. For the IW function, the simplest form is assumed, namely 

t(y) = l-p 2 (y-l) . (13) 

For the heavy-light form factors, LINSR uses those computed with light-cone 
QCD sum rules p 



• NRSX : the second model is the one defined in ref. [19] and makes use of the 
functions f3i(y) and ji(y) calculated in refs. |L(| and 17], respectively. The IW 
function is obtained using a relativistic oscillator model, which gives 

fM = JLlM-^}> d4) 



3 The computation of the 1/m corrections is model dependent, relying on the evaluation of a set of hadronic 
matrix elements of higher dimensional operators in the HQET. 

6 



where A = 2p 2 — 1. Concerning the heavy-light form factors, NRSX improves the 
old WSB model [ 20 1 by implementing, for q 2 ~ g^ m , the expected heavy quark 
scaling laws, see ref. |T^] for details. 

We end this section by spending a word of comment about the parameter p 2 entering 
eqs. ( |l3| ) and (|l4|). In both models, p 2 is defined as the slope of the IW function at the 
zero-recoil point {i.e. y = vb ■ v = 1), namely 



d_ 

dy 



(15) 



y=l 



The value of p 2 is related to the semileptonic differential rate for B — > D*lv, 



dT(B -> D*lv) 
dy 



G 



48vr 



F M D ,(M B -M D *) 2 x/r 



i + 



4y 1 — 2yr + r 

y + l (i-i 



2l 



M^y) 2 



(16) 



where r = Mb* /Mb and ^(y) is an effective semileptonic form factor. The latter is a 
calculable function of the £j(y)s. To make contact with experiments, one defines the 
slope 



P 2 - 



1 d 



Hv) 



y=l 



(17) 



HI) dy' 

which can be extracted from the measurement of the semileptonic differential rate. 

The relation between the experimental slope p 2 and the theoretical parameter p 2 
depends on the model used for the calculation of eq. (|l6|) . We found 



LINSR 
NRSX 



P 2 



p 2 - 0.13, 

p 2 + 0.084, (18) 
2 (0.14 ± 0.02) + 0(l/m) and 



P 



to be compared with the results of ref. [21], p 2 
p 2 = p 2 ± 0.2 respectively. 

In the next section, we will study the dependence of our results on the physical 
slope p 2 , rather than p 2 . 



3 Fitting Cabibbo- allowed decay modes 

This section contains the results of our phenomenological analysis of Cabibbo- 
allowed B decays, focused on the role of heavy-heavy form factors. We proceed as 
follows: 

• we show that the best fit to the BR{B — > Dir) and BR(B — > KJ/^f) is obtained 
for different values of p 2 , depending on the model used for computing heavy-heavy 
form factors; 

• we use the ratios R-w{B — > Dir) introduced in ref. [24], see below, and show 
that, by using them instead of the BR(B — ► Dir), it is possible to fit the FP 
(almost) independently of p 2 . This method gives our best determination of the 
FP, free from theoretical uncertainties coming from the assumptions made for the 
heavy- heavy form factors; 
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Channel 


Si? xlO s 




*B d - 


-> 7T+D- 


300 ± 40 




*B d - 


-> ir + D*- 


276 ± 21 




*B d - 


-> p + D- 


790 ± 140 




*B d - 


-> 


670 ± 330 




B d ^ 


£>+£T 


800 ± 300 


Type I 


B d ^ 


D+D*- 


960 ± 340 


oc \a± 2 


B+ - 


D+D° 


1300 ± 400 




B+ - 


D+D*° 


1200 ± 500 




B d ^ 


D*+D~ 


1000 ± 500 




B d ^ 


D*+D*- 


2000 ± 700 




B+ - 


D*+D° 


900 ± 400 




B+ - 


D* + D*° 


2700 ± 1000 




B d ^ 




89 ± 12 


Type II 


B d ^ 


K*°J/ip 


135 ± 18 


oc a 2 2 


B+ - 


K + J/ip 


99 ± 10 




B+ - 




147 ± 27 




*B+ 




530 ± 50 


Type III 


*B+ 


-» vr+D*° 


460 ± 40 


oc \x\a\ + X2C12 2 


*B+ 




1340 ± 180 




*5+ 


-» 


1550 ± 310 



q 2 


j^BR{B -» D*l+v)(GeV-' 2 ) 


-£pBR(B -► Dl+u)(GeV-' 2 ) 


5 


(0.237 ± 0.026) x 10~ 2 


(0.35 ± 0.06) x 10" 2 


2 

m P 


(0.250 ± 0.030) x 10~ 2 


(0.33 ± 0.06) x 10~ 2 


m D s 


(0.483 ± 0.033) x 10~ 2 




m DZ 


(0.507 ± 0.035) x 10~ 2 





Table 3: Experimental branching ratios \2^, of decay modes to be used in the fit of 
the parameters a, £ and 5%. The classification of non-leptonic channels according to their 
dependence on a± and a 2 is also shown. The channels marked with * have been used for the 
fit of the ratios defined in eq. ($W. 
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• using the (/^-independent) FP, we perform a fit to the BR(B — > Dir) in order to 
extract a preferred range for the value of p 2 ; the results are model dependent, in 
agreement with our first finding; 

• we show that the different ranges of p 2 actually correspond to the same values of 
the relevant form factors at q 2 = M 2 ~ 0. Using the HQET, the latters can be 
determined by factorization applied to B — > Dir decays and may be compared 
with direct measurements from semileptonic decays. 

The relevant decay modes which we use in the fits are listed in tab. || It is well 
known that, in the factorization approximation, only two combinations of DE and CE 
appear in the amplitudes of these decays. This feature is taken into account by the 
parameters a\ and introduced by BSW || . The relation between a\ and and our 
parameters is given by 

ax = a (C 2 + £e^Ci) , fl2 = a (Ci + ^C 2 ) (19) 

where C\ and C2 are the Wilson coefficients of the operators defined in eq. (|2|) i. In 
tab. [31 the decay modes are organized according to the standard classification in three 
classes. Amplitudes of Type I, II and III decay modes are proportional to |ai|, \aa\ 
and \x\a\ A- X2d2\ respectively, where Xi are generic, process-dependent coefficients. 

Given the structure of the amplitudes, we have to fit decay modes of all the three 
classes in order to fully determine the FP. Note that the three classes have a different 
dependence on the form-factors. While Type-II decays always involve a heavy-light 
transition, heavy-heavy form factors enter Type-I modes only. The latter is a general 
feature, since Type-I transitions are always driven by charged currents and are therefore 
proportional to a%. In general Type-III modes involve transitions of both sorts. In our 
parameterization, Type-I modes essentially fix a, while both Type II and III are needed 
to constrain £ and 5%. 

Contrary to the common wisdom, the assumptions made on the momentum de- 
pendence of the heavy-heavy form factors introduce large uncertainties in the deter- 
mination of the FP from this fit. This is more easily shown fitting only the effective 
number of colour as done in the old literature 13. In our parametrization this corre- 
sponds to assume a = 1 and 5% = and to fit £ only. The result of this fit is shown 
in fig. [I], where minimum values of X 2 /dof from the fit of £ to the BR(B —> Dir) and 
to the BR(B — * KJ/^) are plotted as function of p 2 . It is apparent that the best 
fit is obtained for quite different values of p 2 , and corresponds to different values of £, 
depending on the model used to compute the heavy-heavy form factors. Consequently, 
in general, the FP fitted using BRs at fixed p 2 , as usually done in the literature, suffer 
from a large theoretical error, which was previously hidden in the choice of a specific 
model when fitting the data. The second important remark is that a comparison of the 
value of p 2 , the "physical" slope measured in semileptonic decays, with that extracted 
from non-leptonic decays is not a good test of factorization since, in the latter case, 
the result is model dependent. 

4 Notice that our operator basis differs from the one of ref. || by the trivial exchange Qi^ *-* Q2,i- 
5 Had we fitted also a, the resulting minimum x 2 /dof would have been almost independent of p 2 , since 

a easily compensates the variation of the form factors with p 2 . Still, the fitted value of a would have been 

strongly p 2 dependent. 
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Figure 1: p 2 dependence of the minimum x 2 /dof from the fit of £ to BR(B — > Dtt) and 
BR(B — > KJ/ty) assuming a = 1 and 5^ = 0. Both LINSR and NRSX results are shown. 



To circumvent this problem in the determination of the FP, instead of the Type-I 
and TypeTII BRs, we fit the ratios 



£ B R(B OWfc)l s3 ,„;7 (20) 

where M = tt, D s , . . . is the emitted meson. The advantage of using eq. ( p0| ) is that in 
these ratios the heavy-heavy form factor dependence drops out completely for Type-I 
decays and is strongly reduced for Type III. In practice, we used the ratios correspond- 
ing to the non-leptonic decays marked with * in table ^. In the fit, besides the ratios 
Rm for Type-I and Type-Ill, we also use all the BRs of the Type-II decays. 

The results of the fit are shown in tab. for NRSX and LINSR. We do not include 

t } (*) 

the D ( *'Ds modes for two reasons: on the one hand, their contribution to the total 
X 2 is suppressed by the large experimental errors in the measured BRs; on the other, 
we want to use them to extract the decay constants fjj s and fo* ■ 

For both choices of form factors, we give the results of two different fits: the first 
includes all types of decays and determines the FP a, £ and 5^. It retains, however, 
a small residual dependence on heavy-heavy form factors, i.e. on p 2 . The second is a 
fit to Type-I and -II channels only, which is totally independent of p 2 . The results are 
quite close. Note that the second fit only involves two combinations of the three FP. 
As a consequence we have to fix one parameter in order to extract the other two: we 
choose to put 8^ = 0, quite consistently with what has been found with the first fit. 
As a consistency check, we have also verified that different values of 5^ do not change 
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LINSR 


/S 2 = 0.80 


/) 2 = 0.90 


-2 i nn 
/r = 1.00 


— > 2>7r) 


X /dof 


1.36 


1.36 


1.36 


Type I 


a 


i no 
1.02 


1.02 


i no 
1.02 


+ 




0.45 


0.45 


0.45 


2522(23 -» fTJ/VO 




0.00 


0.00 


0.00 


2^(5 -» Dtt) 




1.40 


1.39 


1.42 


Type I+III 


a 


1.05 


1.04 


1.03 


+ 




0.44 


0.44 


0.44 


BR(B -» ifjyv) 




0.00 


0.00 


-0.26 



NRSX 


/5 2 = 1.25 


p 2 = 1.35 


p 2 = 1.45 


22* (J3 Dir) 




0.39 


0.39 


0.39 


Type I 




1.01 


1.01 


1.01 


+ 




0.38 


0.38 


0.38 


2522(23 2iTJ/^) 




0.00 


0.00 


0.00 


2^(23 -» L>tt) 




0.71 


0.70 


0.68 


Type I+III 


a 


1.04 


1.04 


1.03 


+ 




0.38 


0.38 


0.38 


BR(B 2TJ/V0 




0.00 


0.00 


0.00 



Table 4: Results of the fit using Dn semileptonic ratios and BR(B — > KJ/^) for three 
different values of the slope p 2 and the two form-factor models, LINSR and NRSX, described 
in the text. For each model, two fits have been performed, the difference being the inclusion 
of Type-Ill Dtt modes, which introduce some p 2 dependence. 

appreciably the results i. In tab. |i| we show the fitted values of the FP for several 
choices of p 2 , for both NRSX and LINSR. As mentioned above, the results turn out to 
be, within a given model, independent of p 2 . 

Having fitted the FP in a /5 2 -independent way, we now use the BR(B — > Dir) to 
extract from the data a preferred range of p 2 . Notice that 5^ is not a critical parameter, 
since the results of the fits are not very sensitive to its value, and that the values of 
a and p 2 are trivially correlated, because the amplitudes only depend on the product 
of a with the effective form factors at q 2 ~ 0. Therefore we choose to perform a 
two-parameter fit of £ and p 2 using the BR(B — ► Dir), at fixed values of 6% and a, as 
extracted from the previous fit of tab. |j. In this way, we can study the correlations 
in the (£,-p 2 ) plane and check the consistency of the determination of £ using different 
fitting procedures. 

Figure |2| shows the contour plots of x 2 /dof in the {£,,p 2 ) plane for NRSX and 
LINSR. The fitted value of £ is consistent with tab. || and the preferred p 2 is larger 
using NRSX than LINSR. Moreover, the LINSR BRs are steeper functions of p 2 , 
consistently with fig. |l[ This observation justifies the choice of the set of values of p 2 
used in tab. ||. 

It is not surprising that the fit to the BR{B — > Dir) gives values of p 2 which are 

6 The fit is not very sensitive to 8%, thus it cannot fix this parameter very precisely, see the final results 
in©. 
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1 . 8 




Figure 2: Contour plots ofx 2 /dof in the plane (£,p 2 ), obtained from a two-parameter fit to 
the non-leptonic B — > Z?7r decays. The other parameters (a and 8%) have been fixed by the 
fit of tab. | in a p 2 -independent way. The solid curves refer to LINSR, the dashed ones to 
NRSX. The different contours correspond to x 2 /dof = (1.5,2,3) x Xmin/dof. 



model dependent. The fit only fixes the values of the relevant heavy-heavy form factors 
fi = f+,A ,... at q 2 = Ml ~ These form factors can be expressed in terms of 
the £jS at q 2 ~ which, in turn, are related to the IW function by heavy quark 
symmetry, see eq. (|l2[). The relation between the fitted form factors at q 2 = M 2 . 
and the values of £i(l), which are fixed by the HQET, depends on the functional form 
adopted for the IW function £(y) and on the value of p 2 . Thus different values of p 2 are 
obtained by fitting the data with different models. In particular, we find that the main 
difference between NRSX and LINSR relies on the choice of the IW function, eqs. ( |i~3| ) 
and (14), rather than in the inclusion of 1/m corrections. Plotting the results of the fit 
in the planes (£, /j), through eqs. (p|)-(|l2|), we obtain almost the same determination 
of the fi(Ml) with both NRSX and LINSR, as shown in figs. | and g. Although we 
have considered only two models in the present study, we believe that this result is 
quite general. 

We stress again that constraints on the heavy-heavy form factors can only be ob- 
tained by combining the results of two independent fits: the first which fixes the FP 
using the ratios R v , that are essentially independent on the model used to calculate 



7 Type-III modes actually depend also on heavy-light form factors, which however appear in color sup- 
pressed contributions to the total amplitude. 
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Figure 3: Determination of various B — ► D7r /orm factors at q 2 = M 2 as functions of £. 
These plots are obtained from fig. using the relations connecting heavy-heavy form factors 
to the Isgur-Wise function, eqs. LINSR (NRSX) form factors are shown in the 

left (right) column. The different contours correspond to \ 2 /dof = (1.5,2,3) x Xmin/dof- 



the form factors, and the second which fixes the form factors using the BR(B — > Dn) 
at fixed values of the FP. 

The comparison of the heavy-heavy form factors directly measured in semileptonic 
decays at q 2 = M 2 ~ with the results in figs. ||| and ^, and table ||, is a real test of 
generalized factorization in B — >• Dtt decays, independently on the choice of the IW 
function and of the value of p 2 . This checks the assumptions we made for computing 
hadronic matrix elements as described in the previous sections. Since we use the 
HQET relations coming from eq. ([l2|), we are left with only one independent form 
factor, namely the IW function. Therefore the only form factor directly measured at 
q 2 = 0, /+(0), already allows a test of our approach. Its value, /+(0) = 0.66 ± 0.06 ± 
0.04 Q, showed as a band in the upper plots of fig. ||, agrees well with the result 
of the fit. The extraction of the other form factors at q 2 = from the CLEO data 
is under way [25] and would allow a check of the HQET relations among B — > D 
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Figure 4: The same plots as in fig. ^for A , Ax and A 2 . 
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form factors near the maximum recoil point. Notice that, at least in principle, the 
fitting procedure described in this section could be used to extract independently the 
values of all the form factors at q 2 = M 2 , just including them among the FP. In this 
case, the comparison of each form factor with the measurements from semileptonic 
decays would be a test of generalized factorization, independent of the HQET relations 
eq. (p^). Unfortunately the present accuracy of the data does not allow a separate 
determination of the different form factors. 

From the discussion above, we conclude that the HQET-inspired parameterization 
of the heavy- heavy form factors in terms of their value at q^ax an d of the slope p 2 , 
which is commonly adopted by experimental collaborations and successfully applied to 
semileptonic decays, is not the most appropriate choice for the factorization analysis 
of non-leptonic decays. 

With the results for the FP given in (|2l|), we can predict BRs of yet-unmeasured 
decay channels, having one D and one light meson in the final state. We list our 
predictions in tab. ^, where the ranges in square brackets give an estimate of the 
theoretical uncertainties. They were found by allowing values of x 2 /dof up to three 
times larger than the minimum. Flavor SU(3) symmetry justifies the use of parameters 
obtained from B —> Dir and B — ► KJ decays to the decays listed in table Large 
flavour effects are unlikely, since the factorized amplitudes already account for some 
SU (3) breaking. 

Finally, we summarize the result of our fits of the FP by quoting their best values 
and ranges of variation, obtained by allowing values of x 2 /dof up to three times larger 
than the minimum. The comparison between the two different models, NRSX and 
LINSR, gives us an estimate of the theoretical uncertainties due to the form-factor 
model dependence. As discussed before, these FP parameters are those obtained using 
the coefficients functions computed at the NLO in the MS scheme with \i = 5 GeV. 
We obtain 





LINSR 


NRSX 


a 


1.04 [0.9-1.2] 
0.44 [0.2-0.5] 
0.0 [-1.5-1.5] 


1.04 [1.0-1.1] 
0.38 [0.2-0.4] 
0.00 [-1.0-1.0] 


l a 2 


1.04 [0.9-1.3] 
0.31 [0.0-0.7] 


1.05 [1.0-1.2] 
0.25 [0.0-0.4] 


x A 

dof 


1.4 [1.4-4.2] 


0.7 [0.7-2.1] 




0.91 [0.8-1.0] 


1.34 [1.1-1.5] 



(21) 



For the sake of comparison with previous literature, we have also shown the values of 
|ai| and \a?\, computed using eq. (|i~9|). It is worth noticing that exact factorization, 
namely a = 1, £ = 1/3 and 5^ = 0, would give values of x 2 /dof 3-4 times larger than 
the fits which use the generalized factorization, for both the models considered here. 

4 Decay constants from DD S decays 

In this section, we extract the meson decay constants fo s and fo* and compare the 
results with available measurements and lattice results. 

We consider the semi- leptonic ratios R Ds ( B -» and R D *{B -» D^D*+), 
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introduced in the previous section, and define the non-leptonic ratios 

D ^ s ' BR{B -> D(*)p+) K ' 

Up to color-suppressed terms, the factorized amplitudes of the decay modes consid- 
ered here are proportional to /d s //tt or to /d*/ f p . Whereas the ratios Rm of eq. (p0|) 
are defined in such a way that the main form-factor dependence drops out, in the 
non-leptonic ratios of eqs. (^), the form factors appearing in the numerator and de- 
nominator are evaluated at different q 2 and do not cancel out. In this case, however, 
it is the dependence on the FP that tends to cancel, as long as penguin contributions 
are neglected. The non-leptonic ratios above are exactly independent of £ and 5^ only 
if charged B + decays are not considered, as done in ref. @. In our case, we prefer 
to double the number of channels in the fit, by including charged B + decays, at the 
cost of introducing a small dependence on FP and on the Z)W decay constants, both 
appearing in color suppressed terms in the decay amplitudes. We take fjj = 200 MeV 
and f D * = 220 MeV @. 

In general, all DD S modes suffer from a further theoretical error. This uncertainty 
originates from using the same FP, obtained from the fit of sec. ||, in the calculation of 
the relevant BRs entering the non-leptonic ratios. Since in DD S decays, the emitted 
meson is heavy, one may expect, according to the LEET approach, larger violations 
to the factorization limit. In other words, in this case the FP may significantly differ 
from those fixed by the Dir and KJ/^> modes. This is a further source of theoretical 
uncertainty, which we are not able to estimate at present. 

Using the the semileptonic ratios R Ds (B -> and R D *(B -> D^D*+) and 

the non-leptonic ratios of eqs. (|22j), the form factors determined from the fit to the 
BR(B — > Dtt) and the FP from (^), we have extracted fu B and fjy*. Results are 
collected in tab. H where we have separately shown the uncertainties coming from the 
experimental errors on the BRs and from the errors on the FP. As before, in order to 
estimate this source of theoretical uncertainty, we present results obtained using both 
LINSR and NRSX. 

From tab. |5], we quote 

f Ds = 270 ± 35 MeV, f D * = 260 ± 30 MeV, (23) 

where the errors indicatively account for all the sources of uncertainty. 

The value obtained for fjj a is in reasonable agreement with the data ||], fo s = 250± 
30 MeV, and with the lattice results f Ds = 218±?2 MeV (quenched), 235t^t^ 7 MeV 
(unquenched) |10|| , although within large experimental and theoretical uncertainties. 
Our prediction for fjj* is in good agreement with the quenched lattice determination, 
f D; = 240 ±20 MeV Q. 

Comparing the NRSX results of tab. [| with the analysis of ref. Q, one finds differ- 
ences of the order of 10-15%. Besides our inclusion of the charged decay modes, this 
difference arises because we take into account contributions from penguin operators, 
which were neglected in ref. . These contributions amount up to 20% in some chan- 
nels, in particular to those used to determine fo a - For this reason, it is worth testing 
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LINSR 


NRSX 


MeV 


semileptonic nonleptonic 


semileptonic nonleptonic 


Id s 
fo* 


304 ± 42 ± 47 253 ± 24 ± 35 
277 ± 36 ± 43 250 ± 31 ± 16 


297 ± 41 ± 26 267 ± 25 ± 21 
274 ± 36 ± 24 261 ± 32 ± 9 



Table 5: Decay constants extracted from both semileptonic and non-leptonic ratios, eqs. ( \2Q ) 
and $2i<\). Both LINSR and NRSX results are shown. The first error comes from the exper- 
imental ones on the BRs, while the second is a "theoretical" error obtained by varying the 
FP in a range corresponding to values of x 2 /dof up to three times larger than the minimum 
one. 



Sf. 



D s 



0.2 
0.1 


-0.1 
-0.2 
-0.3 
-0.4 
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0.6 
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Figure 5: Charming penguin effect on the determination of fn s from semileptonic ratios 
using NRSX. The difference 5frj a = fD 3 (VL, ^l)//d s (0, 0) — 1 is plotted as a function of rj L 
at different values of 5l- 



the effect of charming-penguin contractions in the determination of the leptonic decay 
constants. We parameterize the effects of charming penguins as in ref. Q, by using 
two real quantities rji and 5l, denoting the relative size and the phase of the charming- 
penguin amplitudes with respect to the corresponding emission ones. In fig. [|, we plot 
Sfo s = fD s {j]Li ^l)//d s (0, 0) — 1 as a function of rjL for various choices of Sl, using 
NRSX. For wl ~ 0.2-0.3 and 5l ~ 7r, as suggested by B — > Kit decays [(jj, |#/d s | is 
about 20%, larger than the (previously) estimated theoretical error on fjj s . Of course, 
there is no compelling theoretical reason to use parameters extracted from Ktt modes 
in this analysis. This exercise shows, however, that penguin effects are not negligible 
and should be included, at least as further source of theoretical uncertainty, at the 



level of 10%, in addition to the one in eq. (|23|). 



5 Conclusions 



In this paper, we have introduced a parameterization of the hadronic matrix elements, 
that generalizes the one of ref. O and express the amplitudes relevant to the calculation 
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of Cabibbo-allowed non-leptonic B decays in terms of factorized matrix elements and 
three real parameters a, £ and 5^. We have shown the connection of our parameteri- 
zation with the generalized factorization of ref. 0. In order to fix these parameters, 
we have reanalized B — > Dtt and B — » KJ/*$> data. 

We have found that the fit to B — > Dtt decays substantially depends on the model 
describing the Isgur-Wise function £ and on the value of its slope. This dependence has 
been drastically reduced by fitting the ratios in eq. (|20|). We have shown that, once that 
the FP are fixed in this way, a best fit to the non-leptonic BR(B — > Dtt) determines 
the values of heavy-heavy form factors at q 2 = M 2 ~ 0. This provides a costraint 
on the HQET models which are currently used for the heavy-heavy form factors. We 
have shown that, in general, different models require different values of p 2 to reproduce 
the fitted values of the form factors. Consequently, a meaningful test of factorization 
is only provided by the comparison of the values of the form factors extracted from 
non-leptonic decays with those directly measured, at small values of q 2 , in semileptonic 
decays. The only form factor directly measured at q 2 = 0, /+(0) = 0.66±0.06±0.04 ||, 
is in good agreement with our finding, suggesting that the generalized factorization 
works well in the case of B — > Dtt decays. 

Our best determination of the FP can be found in (^l|), where an estimate of the 
theoretical uncertainties is also given. Using these FP, we have also presented a set of 
predictions for the BRs of yet-unmeasured B decays, including Dtt, Dp, DK^ modes, 
see tab. ^. 

Finally, using non-leptonic ratios of eqs. ( p0| ) and (p2|), we have extracted the 
charmed meson decay constants from the BR{B — > DD S ), finding 

f Ds = 270 ± 45 MeV, f D % = 260 ± 40 MeV, (24) 

where errors indicatively account for all sources of uncertainty present in the fit, in- 
cluding charming penguins. 

Acknowledgements 

We warmly thank L. Lellouch for very useful discussions. M. Artuso kindly informed 
us about the current status of CLEO analyses of the B — > D^ form factors. We 
acknowledge helpful comments on the experimental results by K.-C. Yang. M.C., E.F. 
and G.M. thank the CERN TH Division for the hospitality during the completion of 
this work. We acknowledge the partial support of the MURST. 

References 

[1] M. Ciuchini, E. Franco, G. Martinelli, and L. Silvestrini, Nucl. Phys. B501 (1997) 
271; 

M. Ciuchini et al, Nucl. Phys. B512 (1998) 3. 

[2] J.-M. Gerard and J. Weyers, UCL-IPT-97-18 ||hep-ph/9711469| ; 
M. Neubert, Phys. Lett. B424 (1998) 152; 

A.F. Falk, A.L. Kagan, Y. Nir, and A. A. Petrov, Phys. Rev. D57 (1998) 4290; 



18 



D. Atwood and A. Soni, Phys. Rev. D58 (1998) 036005; 

for a review, see e.g. R. Fleischer, Int. J. Mod. Phys. A12 (1997) 2459. 

L. Maiani and M. Testa, Phys. Lett. B245 (1990) 585. 

A.J. Buras, J.-M. Gerard, and R. Ruckl, Nucl.Phys. B268 (1986) 16. 

M.J. Dugan and B. Grinstein, Phys. Lett. B255 (1991) 583. 

M. Bauer, B. Stech, and M. Wirbel, Z. Phys. C34 (1987) 103. 

M. Neubert and B. Stech, CERN-TH/97-99, Hhep-ph/9705292(| , to appear in 



"Heavy Flavours II", edited by A.J. Buras and M. Lindner. 
T. Bergfeld et al (CLEO Coll.), CLEO CONF 96-3. 

A. Nippe, to appear in the proceedings of the 7th International Symposium on 
Heavy Flavor Physics, Santa Barbara, California, July 1997. 

T. Draper, plenary talk at "Lattice '98", 13-18 July 1998, Boulder, Colorado, 
USA. 

D. Becirevic et al, ROME1-1227/98, in preparation. 

J.D. Bjorken, Nucl. Phys. B (Proc. Suppl.) 11 (1989) 325. 

M. Bauer and B. Stech, Phys. Lett. B152 (1985) 380; 
F. Buccella et al, Phys. Rev. D51 (1995) 3478. 

A.J. Buras, Nucl. Phys. B434 (1995) 606; 

A.J. Buras and L. Silvestrini, TUM-HEP-315/98, |hep-ph/980627^| . 

M. Ciuchini et al, Z. Phys. C68 (1995) 239. 

M. Neubert and V. Rieckert, Nucl. Phys. B382 (1992) 97. 

M. Neubert, Nucl. Phys. B371 (1992) 149. 

V.M. Belayev et al, Phys. Rev. D51 (1995) 6177; 

P. Ball, Phys. Rev. D48 (1993) 3190; 

P. Ball and V.M. Braun, Phys. Rev. D54 (1996) 2182. 

M. Neubert, V. Rieckert, B. Stech, and Q.P. Xu, in "Heavy Flavours", edited by 
A.J. Buras and M. Lindner (World Scientific, Singapore, 1992), p. 286. 

M. Wirbel, B. Stech, and M. Bauer, Z. Phys. C29 (1985) 637. 

M. Neubert, Phys. Lett. B338 (1994) 84; 

M. Neubert, proceedings of the 27th ICHEP, Glasgow, Scotland, Jul 20-27, 1994, 
edited by P.J. Bussey and I.G. Knowles. IOP, 1995. 

C. Caso et al (PDG), Eur. Phys. J. C3 (1998) 1. 

T. Browder and K. Homscheid, Prog. Nucl. Part. Phys. 35 (1995) 81. 

J. .D. Bjorken, Nucl. Phys. B (Proc. Suppl.) 11 (1989) 325. 

M. Artuso, private communication. 

J.M. Flynn and C.T. Sachrajda, SHEP 97/20, fliep-rat/971005711 , to appear in 



"Heavy Flavours II", edited by A.J. Buras and M. Lindner; 
R.M. Baxter et al. (UKQCD Coll.), Phys. Rev. D49 (1994) 474. 



19 



